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^ ' Abstract. First preliminary results of the balloon-borne experiment SPHERE-2 on the all- 

nuclei primary cosmic rays (PCR) spectrum and primary composition are presented. The 
primary spectrum in the energy range 10^®-5 • 10^^ eV was reconstructed using characteristics 
^ ' of Vavilov-Cherenkov radiation of extensive air showers (EAS), reflected from a snow surface. 

, Several sources of systematic uncertainties of the spectrum were analysed. A method for 

separation of the primary nuclei' groups based on the lateral distribution function' (LDF) 
\ steepness parameter is presented. Preliminary estimate of the mean light nuclei' fraction /30-150 

at energies 3 • 10^'^-1.5 • lO" eV was performed and yielded /3o-i5o= (21±11) %. 
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1. Introduction 

Despite more than 50 years of the cosmic rays' (CR) energy spectrum measurements, a 
considerable uncertainty of the spectral shape in the ultrahigh energy domain (E > 10^^ eV= 1 
5^ . PeV) still exists (for recent results compilation see [1], fig. 15). The difference between spectra 

reconstructed by the various experimental groups is sometimes quite pronounced, and cannot be 
attributed to the statistical issues only. For example, the CR all-nuclei spectrum measured by 
the GAMMA group [2j exhibits statistically significant (> 4(t) sharp peak at 8-10^^ eV. A similar 
feature exists in the spectrum based on the 2009-2010 run of the Tunka-133 experiment [3]~[1]) 
but for the 2010-2011 run of the same experiment it was not confirmed [3]. Also, the peak is not 
seen by other CR experiments such as Akeno [5], Tibet-Ill [6j and KASCADE-Grande [7]-[8]. 

The situation with the primary CR composition at superhigh energies is even more enigmatic, 
as the consensus between the different measurements, such as [9]-|12j. is virtually absent (see 
fig. 1 of [TS])- More recent review J5| claims that the results on the mean logarithmic mass 
number, < InA >, are still consistent with a trend of growing at energies E >3-4 PeV, where 
the primary CR spectrum exhibits a pronounced break, the so-called knee [15]. But anyway, 
the scatter of < InA > values measured by the different experiments is high, and more detailed 



information than just < InA > is needed for the detailed comparison with the CR acceleration 
and propagation models. 

These examples highlight the importance of the all-nuclei CR spectrum measurement and 
composition studies in various experiments and by the diverse methods. In present work 
we present the first results of the SPHERE-2 experiment. The SPHERE-2 experiment is a 
Cherenkov telescope situated above the snow surface of Lake Baikal and registers reflected 
Vavilov-Cherenkov radiation ("Cherenkov light"). The basic idea of such an experiment was 
proposed in [TB]. 

2. The data 

The SPHERE-2 experiment uses a mosaic of 109 PMTs for recording the light signal reflected 
from the spherical mirror. Observation level is typically 400-900 m above the snow surface. In 
the present work we use the data from the 2011 and 2012 experimental runs. Time sampling 
step is 25 ns for the 2011 run and 12.5 ns for the 2012 run. For more detailed description of the 
SPHERE-2 experiment see tl7j-[18j. 

In the four flights of the 2011 run 230 events consistent with EAS Cherenkov light reflected 
from the snow surface were found. For the 2012 run we use the data taken in the first four flights, 
and the number of events classified as showers is 301. Relatively small number of the showers 
compared to the ground-based arrays is due to a high threshold of the EAS observation, as, in 
average, only 1 out of each ~10^ Cherenkov photons reflected from the snow surface reaches the 
detector's mirror. 

3. A method for the primary spectrum and composition reconstruction 

In order to reconstruct the primary CR spectrum and composition we applied an approach based 
on detector response simulation by means of the full Monte Carlo method. The procedure of 
the all-nuclei spectrum measurement consists of the several stages. 

1. As the first step, the full Monte Carlo simulation of the lateral distribution function (LDF) 
of the EAS Cherenkov light, as well as the time structure of the shower was performed using the 
standard code CORSIKA6.500 [191 (for more details see [20j-|22]). All results presented here 
are based on the QGSJET-I high energy hadronic model |23j and the GHEISHA low energy 
hadronic model |24] . 

2. Then, the geometrical and optical effects smearing the observed signal were accounted 
for. For the 2011 run reconstruction we use the simplified detector response model developed 
in |21j-|22j and written in C-I--I-. In the case of the 2012 run the more advanced model is used, 
which was implemented using the standard code Geant4 [25]. 

3. The simulation of the instrumental acceptance was performed. The result is most sensitive 
to the observation altitude's change, as well as to the variations of the thresholds of the detector's 
channels. 

4. For the sample of the EASs registered in the experiment the primary energy was estimated. 
We use the method of [26j, where the energy estimation was performed together with the axis 
position reconstruction by normalising of the experimental LDFs to the model LDFs with known 
energy. 

5. The all- nuclei CR spectrum was reconstructed using the sample of the showers with 
estimated energy and results of the acceptance calculation. 

The primary composition study in the present work is based on the LDF steepness parameter, 
defined as the ratio of the number of Cherenkov photons in the circle with the radius of 67 m 
to the same number in the ring with the radii of 67 m and 134 m [21]. The parameter allows to 
select some fraction of the light nuclei [2D]-[5T]. In present work protons stand for light nuclei, 
and Iron for heavy ones. 



4. Preliminary results 

The all-nuclei CR spectrum reconstructed using 2011-2012 runs data of the SPHERE-2 
experiment is shown in fig. 1 (red stars). Statistical (red bars) and systematic (red dashed 
lines) uncertainties of the spectrum are shown. The main sources of systematic uncertainty 
of the spectrum are as follows: a) bin-to-bin migration due to the energy estimation error 
(dominant factor at E>20 PeV); b) acceptance estimation inaccuracy (not more than 5 %); c) 
the spectal shape vs. the primary composition dependence (dominates at E<20 PeV). With 
respect to the direct predecessor of the present work [27], we achieved a considerable progress 
of the reconstruction methodology. Also the results of some other experiments are shown in 
fig. la for comparison: Akeno [3] (green circles), KASCADE-Grande [7] (black triangles), and 
Tunka-133 [28j (blue squares). Statistical errors of the Akeno spectrum are comparable with the 
circle's diameter; systematic uncertainties of the KASCADE-Grande result are shown by black 
dashed lines. 

Estimated fraction of light nuclei / vs. lg{E) is shown in fig. 2 (red circles) with statistical 
uncertainties of the / values (red bars). Several independent reconstructions of the experimental 
showers' LDFs, and then the / vs. lg{E) dependence were performed, and induced systematic 
uncertainty of the / value was estimated (black lines). In the first two bins the value of / is 
strongly modified by the threshold effects; an additional uncertainty assotiated with them is 
shown in the picture by blue arrows. The data presented in fig. 2 makes it possible to estimate 
the mean (not intensity- weighted!) fraction of light nuclei in the energy region 30-150 PeV 
/30-150 = 0.21 lb 0.11. The corresponding value of < InA >= 3.20 it 0.45 (comparable with 
values allowed by [14j; see fig. 17). 

All presented results are subject to additional recheck and develompent of the methodology, 
and the detailed study of the systematic uncertainties that might change the results presented 
at figs. 1-2 is underway. 
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Figure 1. (color online). The all- nuclei 
CR spectrum measured in the SPHERE-2 ex- 
periment (red stars; preliminary result) to- 
gether with the statistical and the system- 
atic uncertainties, as well as results of the 
Akeno [5] (green circles), KASCADE-Grande 
[71 (black triangles) and Tunka-133 [28] exper- 
iments (blue squares). 
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Figure 2. (color online). The reconstructed 
fraction of light nuclei / vs. lg{E), as 
well as the statistical (red bars) and the 
systematic (black lines) uncertainties of /. 
Also the additional error of the light-nuclei 
fraction arising from the threshold effects 
(blue arrows) is shown. 



5. Conclusions 

In the present work we performed the first detailed reconstruction of the CR all-nuclei spectrum 
using Cherenkov light reflected from the snow surface, as well as the very first study of the 
primary composition by means of this method. Considerable effort was devoted to establishing 
the database of the full Monte Carlo model showers in order to correctly simulate their 
fluctuations. A dependence of the spectral shape vs. primary composition in the threshold 
region was accounted for. It was found that the overall fraction of light nuclei in the energy 
region 30-150 PeV /30-150 = 0.21 it 0.11. All above-mentioned results are preliminary and are 
subject to the additional recheck. 

To conclude, let us note that the developed method of the all-nuclei spectrum reconstruction 
and composition studies might be useful for analysing of the next generation CR experiments 
capable of observation of Cherenkov light reflected from the natural surfaces, particularly, for 
orbital experiments such as TUS [29j and JEM-EUSO |30j . 
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